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qq ' Abstract 
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o 



Theoretical understanding of experimental results from relativistic 
heavy-ion collisions requires a microscopic approach to the behav- 
ior of QCD n-point functions at finite temperatures, as given by the 
OO ■ hierarchy of Dyson-Schwinger equations, properly generalized within 

the Matsubara formalism. The technical complexity of related finite- 
temperature calculations however mandates modeling. We present a 
■ model where the QCD interaction in the infrared, nonperturbative do- 

main, is modeled by a separable form. Results for the mass spectrum 
f4 of light quark flavors at finite temperature are presented. 

1 Introduction 

> ■ 

^ ■ While the experiments at RHIC [U [2] advanced the empirical knowledge of 

the hot QCD matter dramatically, the understanding of the state of matter 
that has been formed is still lacking. For example, the STAR collaboration's 
assessment [2] of the evidence from RHIC experiments depicts a very in- 
tricate, difficult-to-understand picture of the hot QCD matter. Among the 
issues pointed out as important was the need to clarify the role of quark- 
antiquark (qq) bound states continuing existence above the critical temper- 
ature T c , as well as the role of the chiral phase transition. 
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Both of these issues are consistently treated within the Dyson- Schwinger 
(DS) approach to quark-hadron physics. Dynamical chiral symmetry break- 
ing (DChSB) as the crucial low-energy QCD phenomenon is well-understood 
in the rainbow-ladder approximation (RLA), a symmetry preserving trunca- 
tion of the hierarchy of DS equations. Thanks to this, the QCD low energy 
theorems are fulfilled and the behavior of the chiral condensate and pion 
mass and decay constant are in accord with the Gell-Mann-Oakes-Renner 
relation, i.e. the Goldstone theorem. 

For recent reviews of the DS approach, see, e.g., Refs. [31 E], of which 
the first [3] also reviews the studies of QCD DS equations at finite tempera- 
ture, started in [5] . Unfortunately, the extension of DS calculations to non- 
vanishing temperatures is technically quite difficult. This usage of separable 
model interactions greatly simplifies DS calculations at finite temperatures, 
while yielding equivalent results on a given level of truncation [6], [7]. A recent 
update of this covariant separable approach with application to the scalar a 
meson at finite temperature can be found in [8]. Here, we present results for 
the quark mass spectrum at zero and finite temperature, extending previous 
work by including the strange flavor. 



2 The separable model at zero temperature 

The dressed quark propagator S q (p) is the solution of its DS equation [31 0], 

S.ip)- 1 = < 7-P + ^ + !/ t0 9 2d °> ~ QlMQl* > (!) 

while the qq' meson Bethe-Salpeter (BS) bound-state vertex T qq -i(p, P) is the 
solution of the BS equation (BSE) 

- \(P 2 )T qq ,(p, P) = \J ^g 2 Df u (p - £) llx S q (£+)r gq ,(l, P)S q (£.) lu , (2) 

where D°^l(p — £) is an effective gluon propagator modeling the nonpertur- 
bative QCD effects, m q is the current quark mass, the index q (or q') stands 
for the quark flavor (u, d or s), P is the total momentum, and £± = £ ± P/2. 
The chiral limit is obtained by setting m q = 0. The meson mass is identified 
from A(P 2 = — M 2 ) = 1. Equations CQ) and (T5]) are written in the Euclidean 
space, and in the consistent rainbow-ladder truncation. 
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The simplest separable Ansatz which reproduces in RLA a nonpertur- 
bative solution of (pQ) for any effective gluon propagator in a Feynman-like 
gauge g 2 Df l/ (p - 1) - 6^D(p 2 , i\ p ■ i) is Cj 

D(p 2 ,l 2 ,p.l) =L>o^o(p 2 )^o(0+^i(P 2 )(p-^i(0 • (3) 

This is a rank-2 separable interaction with two strength parameters Di and 
corresponding form factors ^(p 2 ), i = 1,2. The choice for these quantities 
is constrained to the solution of the DSE for the quark propagator ([I]) 

S.ip)- 1 = i 1 ■ pA q (p 2 ) + B q (p 2 ) = Z-^p 2 )^ ■ p + m q (p 2 )} , (4) 

where m q {p 2 ) = B q {p 2 ) / A q {p 2 ) is the dynamical mass function and Z q {p 2 ) = 
A^ip 2 ) the wave function renormalization. Using the separable Ansatz ([3]) 
in ([1]), the gap equations for the quark amplitudes A q {p 2 ) and B q {p 2 ) read 

Once the coefficients a q and b q are obtained by solving the gap equations (J5]) 
and (EI), the only model parameters remaining are rh q and the parameters of 
the gluon propagator, to be fixed by meson phenomenology. 



3 Extension to finite temperature 

The extension of the separable model studies to the finite temperature case, 
T 7^ 0, is systematically accomplished by a transcription of the Euclidean 
quark 4-momentum via p — ► p n = (u n ,p), where u n = (2n + 1)ttT are the 
discrete Matsubara frequencies. In the Matsubara formalism, the number 
of coupled equations represented by ([TJ and (j2J) scales up with the num- 
ber of fermion Matsubara modes included. For studies near and above the 
transition, T > 100 MeV, using only 10 such modes appears adequate. Nev- 
ertheless, the appropriate number can be more than 10 3 if the continuity 
with T = results is to be verified. The effective qq interaction will auto- 
matically decrease with increasing T without the introduction of an explicit 
T-dependence which would require new parameters. 
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The solution of the DS equation for the dressed quark propagator now 
takes the form 



S q 1 {Pn, T) =ij -p A q (p 2 n , T) + i^u n C q {p 2 n , T) + B q (p 2 n , T), (7) 

where p 2 n = u 2 +p 2 and the quark amplitudes B q (p 2 l ,T) = m q + b q (T)J r o(p 2 l ), 
A q {p 2 n ,T) = l + a^T^ipl), and C q {p 2 n ,T) = 1 + c q {T)F x (p 2 n ) are defined 
by the temperature-dependent coefficients a q (T),b q (T), and c q {T) to be de- 
termined from the set of three coupled non-linear equations 

<T) = J ^- 3 ^(p 2 n )p 2 [l + a q (T)^(p 2 n )} d~\plT) , (8) 

c q (T) = / ^Hpl)^ + o q (T)Hp 2 n )] d~\plT) , (9) 

b q (T) = ^TY: J ^MpD [™ q + \(T)T Q (pl)\ d-\plT) . (10) 

The function d q {p 2 n , T) is the denominator of the quark propagator S q (p n , T), 
and is given by 

d q (p 2 n , T) = p 2 A 2 q (p 2 n , T) + u 2 C 2 (p 2 n , T) + B 2 (p 2 n , T). (11) 

The procedure for solving gap equations for a given temperature T is the 
same as in T = case, but one has to control the appropriate number of 
Matsubara modes as mentioned above. 



4 Confinement and Dynamical Chiral Sym- 
metry Breaking 

If there are no poles in the quark propagator S q (p) for real timelike p 2 then 
there is no physical quark mass shell. This entails that quarks cannot prop- 
agate freely, and the description of hadronic processes will not be hindered 
by unphysical quark production thresholds. This sufficient condition is a 
viable possibility for realizing quark confinement [!7J. A nontrivial solution 
for B q (p 2 ) in the chiral limit (m = 0) signals DChSB. There is a connection 
between quark confinement realized as the lack of a quark mass shell and the 
existence of a strong quark mass function in the infrared through DChSB. 
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The propagator is confining if m 2 (p 2 ) ^ —p 2 for real p 2 , where the quark 
mass function is m q {p 2 ) = B q (p 2 )/A g (p 2 ). In the present separable model, 
the strength b q = B q (0), which is generated by solving Eqs. (jSJ) and (jSJ), con- 
trols both confinement and DChSB. At finite temperature, the strength b q (T) 
for the quark mass function will decrease with T, until the denominator (JTTj) 
of the quark propagator can vanish for some timelike momentum, and the 
quark can come on the free mass shell. The connection between deconfine- 
ment and disappearance of DChSB is thus clear in the DS approach. Also 
the present model is therefore expected to have a deconfinement transition at 
or a little before the chiral restoration transition associated with bo(T) — > 0. 
The following simple choice of the separable interaction form factors, 

T I 2\ / 2/a2\ t~ I 2\ 1 + ex P( — Po/^-l) 

Fo[P = exp(-p /A , fi p = — — 2V772T > 

1 + exp((p 2 -pg)/Af) 

is used to obtain numerical solutions which reproduce very well the phe- 
nomenology of the light pseudoscalar mesons and generate an acceptable 
value for the chiral condensate. 

The resulting quark propagator is found to be confining and the infrared 
strength and shape of quark amplitudes A q {p 2 ) and B q {p 2 ) are in quantita- 
tive agreement with the typical DS studies. Gaussian-type form factors are 
used as a minimal way to preserve these properties while realizing that the 
ultraviolet suppression is much stronger than the asymptotic fall off (with 
logarithmic corrections) known from perturbative QCD and numerical stud- 
ies on the lattice [9]. 

5 Results 

Parameters of the model are completely fixed by meson phenomenology cal- 
culated from the model as discussed in [3 [8] . In the nonstrange sector, we 
work in the isosymmetric limit and adopt bare quark masses m u = = 5.5 
MeV and in strange sector we adopt m s = 115 MeV. Then the parameter 
values 

A = 758 MeV, Ai = 961 MeV, p = 600 MeV, (12) 

A)Ao = 219, DiA^ = 40, (13) 

lead, through the gap equation, to a u ^ d = 0.672, b u d = 660 MeV, a s = 0.657 
and b a = 998 MeV i.e., to the dynamically generated momentum- dependent 
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mass functions m q (p 2 ) shown in Fig. [TJ In the limit of high p 2 , they converge 
to m u and m s . However, at low p 2 , the values of m u (p 2 ) are close to the 
typical constituent quark mass scale ~ m p /2 ~ 171^/3 with the maximum 
value at p 2 = 0, m u (0) = 398 MeV. The corresponding value for the strange 
quark is m s (0) = 672 MeV. Thus, Fig. [1] shows that in the domain of 
low and intermediate p 2 < 1 GeV 2 , the dynamically generated quark masses 
have typical constituent quark mass values. Thus, the DS approach provides 
a derivation of the constituent quark model [10J from a more fundamental 
level, with the correct chiral behavior of QCD. 

Obtaining such dynamically generated constituent quark masses, as previ- 
ous experience with the DS approach shows (see, e.g., Refs. such as [31IH[TU]). 
is essential for reproducing the static and other low-energy properties of 
hadrons, including decays. (We would have to turn to less simplified DS 
models for incorporating the correct perturbative behaviors, including that 
of the quark masses. Such models are amply reviewed or used in, e.g., Refs. 
[31 IH QUI EH] , but addressing them is beyond the present scope, where the 
perturbative regime is not important.) 




Figure 1: The p 2 dependence (at T = 0) of the dynamically generated quark 
masses m s (p 2 ),m u (p 2 ) for, respectively, the strange and the (isosymmetric) 
nonstrange case. 

The extension of results to finite temperatures is given in Figs. [21 [3j 
Very important is the temperature dependence of the chiral-limit quantities 
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Bq(0)t and (qq)o(T), whose vanishing with T determines the chiral restora- 
tion temperature T Ch . We find T Ch = 128 MeV in the present model. 
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Figure 2: The temperature dependence of B s (0), B u (0) and B (0), the scalar 
propagator functions at p 2 = 0, for the strange, the nonstrange and the chiral- 
limit cases, respectively. The temperature dependence of the chiral quark- 

1 /3 

antiquark condensate, — (qq) , is also shown (by the lowest curve). Both 
chiral-limit quantities, -Bo(0) and —(qq)^ 3 , vanish at the chiral- symmetry 
restoration temperature Tch = 128 MeV. 

The temperature dependences of the functions giving the vector part of 
the quark propagator, A u ^ s (0)t and C UjS (0)t, are depicted in Fig. @J Their 
difference is a measure of the 0(4) symmetry breaking with the temperature 
T. 

The presented model, when applied in the framework of the Bethe-Salpeter 
approach to mesons as quark- antiquark bound states, produces a very satis- 
factory description of the whole light pseudoscalar nonet, both at zero and 
finite temperatures [12]. 
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Figure 3: The temperature dependence of m s (0),m u (0) and m (0), the dy- 
namically generated quark masses at p 2 = for the strange, the nonstrange 
and the chiral-limit cases, respectively. 
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Figure 4: The violation of 0(4) symmetry with T is exhibited on the example 
of A U , S (0) T and C U , S (0) T . 
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